Extrinsic models for the dielectric response of CaCu3Ti40i2 
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The large, temperature-independent, low-frequency dielectric constant recently observed in single- 
crystal CaCu3Ti40i2 is most plausibly interpreted as arising from spatial inhomogenities of its local 
dielectric response. Probable sources of inhomogeneity are the various domain boundaries endemic 
in such materials: twin, Ca-ordering, and antiphase boundaries. The material in and neighboring 
such boundaries can be insulating or conducting. We construct a decision tree for the resulting 
six possible morphologies, and derive or present expressions for the dielectric constant for models 
of each morphology. We conclude that all six morphologies can yield dielectric behavior consistent 
with observations and suggest further experiments to distinguish among them. 
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I. INTRODUCTION 

The low-frequency hmit eq of the dielectric constant 
e{uj) of CaCu3Ti40i2 (CCTO) single-crystals can be as 
much as 10^ times larger than that expected by extrap- 
olating from electron and infrared phonon contributions, 
approaching 10^ and remaining constant over a broad 
temperature range. ^"'^ With a decrease in temperature 
or increase in frequency, e falls off, displaying Debye-like 
relaxation behavior with an activated relaxation rate.^''^ 
Such unexpected behavior is startling scientifically and 
intriguing technologically and, accordingly, has attracted 
much attention,*"^ including comparative studies involv- 
ing its Cd counterpart, CdCu3Ti40i2 (CdCTO), which 
apparently shows similar though less pronounced behav- 
ior in ceramic samples. ^'^"^ 

The central question is now whether the large dielectric 
response is intrinsic to a perfect crystal of CCTO or ex- 
trinsic in that it originates with defects, inhomogeneities, 
etc., in particular samples. Based on our first-principles 
calculations of the structure, the phonon spectrum, and 
the dynamical effective charges of both CCTO (Ref. 7) 
and CdCTO (Ref. 8) , together with the agreement of our 
results with experiment^'-'^'' and an assessment of the ex- 
isting experimental facts, ^^^^^^ we have argued strongly 
against an intrinsic interpretation, as have others on em- 
pirical grounds.^'^^'^^ At issue now is the nature of the 
extrinsic mechanism, which we address in the present pa- 
per. 

Empirically, the single-crystal samples of CCTO are 
known to be to be highly twinned, "'^ and the transport 
behavior of these domains and their boundaries, as well 
as those of other such domains and boundaries (in both 
single-crystal and ceramic samples), could play a signif- 
icant role in the observed low-frequency giant dielectric 
response, as was noted in earlier work.^'*"'^^ The basic 
idea is that the bulk of the material is either conduct- 
ing or nearly so, and that the conductivity of the en- 
tire sample is only prevented either by a failure of the 
conducting regions to percolate, or else by the presence 



of thin insulating blocking layers at the surfaces or at 
internal domain boundaries. The various morphologies 
associated with these possibilities can all be consistent 
with an enormous enhancement of the static dielectric 
constant and with a Debye-like frequency response, al- 
though there are characteristic details of the dielectric 
response that may help distinguish between them. 

In this article, we identify and thoroughly analyze 
those morphologies consistent with existing experimen- 
tal results for CCTO. We begin in Sec. II by listing six 
morphology classes, conceived from general arguments, 
that are associated with internal and external bound- 
aries probably present in CCTO and may give rise to 
the observed dielectric response. Internal domain bound- 
aries present in CCTO could be conducting or noncon- 
ducting, and if the former is the case, the bulk must 
be insulating for a large dielectric response to be possi- 
ble; alternatively if the latter is the case, the bulk must 
be conducting. Of course, given either instance, intrin- 
sic or electrode-induced insulating (or blocking) layers 
must be present at the surface if the internal conduct- 
ing regions percolate through the sample, since a finite 
dielectric response is observed. And, if they do not per- 
colate, blocking layers may or may not be present. Af- 
ter sorting out the possible morphological scenarios, the 
dielectric response of blocked morphologies is discussed 
in Sec. III. We then provide useful exact bounds on eo 
in Section IV, and afterward analyze the two unblocked 
morphologies in detail, completing our analysis. The di- 
electric properties of a solid exhibiting bulk conductivity 
with insulating boundaries and no blocking are derived 
in Section V; the converse morphology is treated in Sec- 
tion VI. We conclude, in Section VII, that any of these 
six morphologies would yield dielectric responses consis- 
tent with observations,^"'^ and that existing experiments 
do not distinguish them. Accordingly we propose various 
experiments, including the use of nanoscale probes to ex- 
amine the local conductivity, to discriminate among the 
various possibilities presented here. 

We do not consider here the possibility that the ex- 
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trinsic behavior arises from point or line defects ran- 
domly but homogeneously distributed throughout the 
bulk. One such possibility would be that oxygen va- 
cancies, common to perovskite-like oxides and which we 
denote Vq, are responsible. This has already been pro- 
posed by Ramirez et al.,^'^ and they point out that a 
Vo concentration of order 10^^ would be required. If 
the activation energies of (Vq)^ and (Vo)^^ i-c, singly- 
and doubly- ionized vacancies arc comparable to those in 
other perovskites^'^^ and if the mobilities of electrons 
in the conduction band are similarly comparable/^"^^ a 
concentration of Vq of order 10~^ would give rise to a dc 
bulk resisitivity many orders of magnitude smaller than 
that observed or inferred for these materials. Accordingly 
we have confined our attention to two-dimensional de- 
fects. 



II. SIX MORPHOLOGIES 

Large quasistatic dielectric constants co can arise in 
macroscopic insulators containing conducting regions ap- 
proaching a percolation threshold. While this may occur 
in principle for a wide range of morphologies, we con- 
sider here those six which we perceive as relevant to the 
good single-crystal CCTO samples possessing the largest 
values of eg.^ Each morphology can be decomposed into 
two or three broadly-defined regions: (1) domains, (2) 
their associated boundaries, and (3) blocking layers (if 
any) parallel to the electrodes. Several domain types 
(and therefore boundaries) are possible. In addition to 
twin boundaries, a necessary by-product of the observed 
twinning, ^'^'^"^ two other categories may exist, both as- 
sociated with variant chemical ordering. The first arises 
from the fact that the Ca sublattice may occupy four 
possible A sites, thus making possible four differently- 
ordered domains and their associated domain boundaries. 
The second includes antiphase domain boundaries of A- 
A or B-B type separating domains in which the A atoms 
have changed registry with the B atoms within the per- 
ovskite structure. Thus, if spatial heterogeneity in the 
electrical response occurs in CCTO single crystals, it 
may be plausibly and economically associated with an 
intersecting set of locally planar twin, antiphase, and/or 
compositional-ordering domain boundaries. 

While each of the boundary regions contain planar 
defects, their composition may differ from that of the 
bulk through the Gibbs adsorption phenomenon; it fol- 
lows that their effective thickness may be substantial in 
the present context. We can thus distinguish two possi- 
bilities. The domains themselves can be conducting or 
insulating in the bulk. If the domains are conducting, 
their boundaries must be blocking, and vice versa. Im- 
purities or oxygen deficiencies could give rise to the bulk 
conductivity. Alternatively, if the impurities or vacancies 
segregate at or near the boundaries, the bulk would be 
insulating and the boundaries conducting. 
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FIG. 1. Morphological decision tree showing how the pos- 
sible morphological elements lead to six distinct morphologies 
which arc nonpcrcolating (NP), percolating (P), conducting 
bulk (CB), conducting interfaces (CI), unblocked (UB), or 
surface blocked (SB). The numbers in parentheses correspond 
to those in Table I. 

Conductivity may also arise at a boundary through the 
formation of bands of interface states, which are localized 
at the boundaries and extended along them, and which 
possess a fundamental gap signficantly smaller than the 
intrinsic bulk gap. This would allow the formation of 
shallow traps associated with the boundary imperfec- 
tions. Conductivity could originate with thermally ac- 
tivated traps or with ordinary thermal excitation across 
the intcrfacial energy gap. For example, consider Cu- 
Cu antiphase boundaries. CCTO is a Mott-Hubbard 
insulator with an intrinsic gap of at least 1.5 eV (Ref. 
3, 7). The highest valence band and lowest conduction 
band states are primarily d-orbitals hybridized with the 
p-orbitals of the four nearest-neighbor oxygens in each 
Cu04 plaquctte.^ The coupling between neighboring pla- 
quettes is weak, as evidenced by the small valence and 
conduction bandwidths. Within a tight-binding model, 
the effective Hubbard U is thus substantially larger than 
the effective number of neighbors z times the effective 
electron transfer integral t {zt is then a measure of the 
bandwidth within this model). However, locally, at a Cu- 
Cu antiphase boundary, z (and therefore the bandwidth) 
increases, U / zt decreases, the Hubbard gap decreases, 
and the possibility of boundary conductivity arises. An 
excess Vq concentration could also give rise to boundary 
conductivity as discussed further in Sec. VII. 

We conclude that in the present state of our knowledge 
of CCTO, we must suppose that the conducting bulk 
(CB) and the conducting interface (CI) morphologies are 
equally plausible. The conducting regions can be below 
the percolation threshold, i.e., nonpcrcolating (NP), or 
above the percolating threshold (P). If the latter (P), 
conduction must be blocked by a surface barrier layer 
(SB). The barrier layer can be intrinsic to the surface or 
be associated with the electrode-CCTO interface. In the 
nonpercolating case, there may or may not be a blocking 
layer. 

These possibilities can be represented by a morpho- 
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TABLE I. Six possible morphologies that could produce 
the large dielectric response. 



and 



MORPHOLOGY 



ACRONYM 



T nonpcrcolating, conducting bulk, unblocked NP/CB/UB 

2 nonpcrcolating, conducting bulk, surface blocked NP/CB/SB 

3 nonperc, conducting interfaces, unblocked NP/CI/UB 

4 nonperc. conducting interfaces, siirface blocked NP/CI/SB 

5 percolating, conducting bulk, surface blocked P/CB/SB 

6 percolating, conducting interfaces, surface blocked P/CI/SB 



logical decision tree, illustrated in Fig. 1. Six distinct 
morphologies result, which are listed in the Table I along 
with the acronyms by which we identify them in the re- 
mainder of the paper. 



III. BLOCKED MORPHOLOGIES 

We suppose that the CCTO samples have plane- 
parallel boundaries normal to their smallest dimension, 
forming good planar electrode-sample interfaces. The ef- 
fective complex dielectric constant e* for all blocked cases 
is then simply 



1 



(1) 



^M 



where /l is the volume fraction of the two blocking lay- 
ers which is equal to the sum of their relative thick- 
nesses, EL their dielectric constant (assumed to be en- 
tirely real), and the complex macroscopic dielectric 
constant of the material. We now focus on two pos- 
sible sample morphologies, both containing conducting 
regions, which lead to a large Dcbye-like low- frequency 
dielectric response; in one the conducting regions perco- 
late, and in the other they do not. 



A. Dielectric response of morphologies with 
percolating conducting regions 

For morphologies in which the bulk material con- 
tains percolating (P) conducting regions, the frequency- 
dependent dielectric constant can be expressed as 



Em = eM(w) + 4ma{u>)/u), 



(2) 



with and a{co) real. Inserting (2) into (1) yields the 
Debye form 



e = Ecx) + 



eo 
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where 
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(4) 
(5) 



Coo 47rc7(w) ' 



(6) 



In (5), eM(oo) is a representative value of eM('^) for iv 
such that ujt{oj) 1. 

It follows from the above that if a percolating mor- 
phology actually occurs, it must be well beyond the per- 
colation threshold. This results from the fact that cm (0) 
diverges at the percolation threshold, so that sufficiently 
close to the threshold one expects no significant differ- 
ence between cq and Cqo; then euioo) > eo, in direct 
contradiction to the observations (eoo -C cq). 

Neglecting e^o, the asymptotic (high-frequency) behav- 
ior of the real part of e*, denoted here e', is 



e ~ 



1 [ATTa{u 



(7) 



thus falling off with uj more slowly than the pure Debye 
case of w"^. Such behavior is also consistent with the ex- 
perimental results of Ref. 3, and thus the two percolating 
blocked cases (P/SB, entries 5 and 6 of Table I), cannot 
be eliminated on the basis of their frequency dependence. 

Using Eq. (4) with a value of el ~ 10^, an eo of order 
10^ can be achieved if we take the blocking layer volume 
fraction as /l ~ 10"'^. As typical samples have thick- 
nesses of the order of millimeters, the implied widths of 
the blocking layers are of the order of microns. It is un- 
likely, though not impossible, that barrier widths are ac- 
tually that large. Moreover, values of eo as low as 5x10^ 
have been observed, implying barrier widths of the or- 
der of tens of microns, which would be improbably large. 
And finally, conducting samples have been fabricated, 
which is inconsistent with the existence of such an elec- 
trode barrier and mitigates against the P/CB/SB and 
P/CI/SB cases. 

In the event that a percolating morphology is demon- 
strated to be the most relevant to CCTO after all, a 
simple experiment would distinguish between intrinsic 
surface and electrode-induced blocking layers. Suppose 
an insulating film of relative thickness fi/2 and dielec- 
tric constant ei is fabricated or deposited between each 
electrode and the sample. The low frequency effective 
dielectric constant of the composite (blocking layers plus 
insulating film), e°g, is given by 



/i _l_ /l 



(8) 



An intrinsic surface blocking layer should be affected but 
not eliminated by the additional insulator. On the other 
hand, the effects of an electrode blocking layer (EB) 
should certainly be eliminated. Thus, measuring l/e^g 
for values of /i greater than the /l for the blocking layer 
and then extrapolating the resulting l/ej^^' vs. fi towards 
zero /i should yield a vanishing intercept for the EB case. 
A similar procedure would yield a finite intercept for the 
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intrinsic case, according to (8). This experiment can be 
carried out by depositing an electrode directly on one face 
of a thin CCTO sample. On the other face, a wedge of 
insulator is deposited followed by the deposition of a set 
of separated electrode strips parallel to the wedge bound- 
ary. A set of values of fi in (10) embracing zero is thus 
achieved with a single sample. 



B. Dielectric response of morphologies with 
nonpercolating conducting regions 

For sample morphologies with nonpercolating conduct- 
ing regions (NP), we use for ej^ the generic form 



IV. EXACT BOUNDS ON THE STATIC 
DIELECTRIC CONSTANT 

Before discussing the nonpercolating imblocked cases 
(NP/CB/UB and NP/CI/UB) and, equivalently, the ori- 
gin of expression (9) for in the nonpercolating blocked 
cases of Sec. TUB, it is convenient to have exact bounds 
on the total static dielectric constant eq. 

We make the simplification that the local dielectric and 
conductivity tensors are isotropic. Suppose an inhomo- 
geneous material has a set {ei} of different local static 
dielectric constant with a corresponding set of volume 
fractions {fi}- The static macroscopic dielectric constant 
eo then lies between the exact Wiener bounds^^ 



1 — iuiTuii^) 



(9) 



The origin of this form will be discussed in more detail for 
the first four cases of Table I in Sees. V and VI. Inserting 
(9) into (1) again yields the Debye form (3) but now with 

/l , 

EL 



en = 



(10) 



and 



eo e 



(11) 



(12) 



For a value of eo close to the observed value of approx- 
imately 10^, it follows from Eq. (10) that fh/^h and 
(1 — /L)/e^ must be each less than 10~^. Since cl should 
be of order 10^, the value of Cb, it follows that /l should 

be 10^"^ or less and e^^ 10^ or larger. For a value of too 
of approximately 10^, it then follows from Eq. (11) that 
eM = e(x)=10^. These two results imply, with (12), that 
t(w) > tm{uj)- 

Interposing insulating films, as in the percolating mor- 
phology in Sec. A above, should not affect the intrinsic 
surface barrier. This implies that 



ei 



1 



(13) 



Cgff CI CM c 

as /i ^ 0. In the electrode-induced case the barrier at 
the electrode is eliminated and /l — > so that 

1 /t l-/i 1 _ 1 



(14) 



In (13) and (14) the arrows signify an extrapolation to 
/i = which is carried out as described in Sec. A above. 

Thus, if the extrapolated value of ejjg- exceeds Cq, one 
learns from this experiment that the sample is blocked 
by an electrode- induced layer and that exceeds eo- 
If e[!jj- does not exceed eo when fi 0, the experiment 
does not distinguish between the two possible sources of 
surface blocking. 



E/4 



H -1 



< e < 



(15) 



the lower bound an equality for all inhomogeneities 
bounded by planes parallel to the electrode planes (con- 
stant electric displacement D) and the upper bound an 
equality for all inhomogeneities bounded by planes per- 
pendicular to the electrode planes (constant electric field 
E). When applied to the imblocked morphologies, these 
bounds are particularly illuminating, as we show below. 



V. UNBLOCKED MORPHOLOGIES WITH 
INSULATING INTERNAL BOUNDARIES 

In the absence of blocking layers, we suppose that there 
is a bulk volume fraction /b with finite conductivity fjb 
and a nonconducting internal interface (boundary) vol- 
ume fraction = 1 — /b with a finite dielectric constant 
ei. Since eb ~ Aniah/co — > oo as w — »■ 0, Eq. (15) simpli- 
fies to 



ei//i < e < 00 . 



(16) 



The \ipper bound pertains to systems above the perco- 
lation threshold. Most interesting, however, is the lower 
bound; a small volume fraction fi ~ 10~^ associated with 
the interfaces and interface dielectric constant of ej com- 
parable to eb together imply a lower bound comparable 
to the observed eo. 

If all internal boundaries are parallel to the electrode 
planes, the frequency-dependent dielectric constant be- 
comes 



ep - eb 
1 — iiVT{u}) 



eo = efc/ fi , 



t{uj) = (1 - /i)eo/47r£76(a;) 



(17) 



(18) 



(19) 
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with /i the volume fraction of the internal blocking lay- 
ers; with the dielectric constant of the blocking layers as- 
sumed equal to eb! and with (Tb{Lo) the bulk conductivity 
within the domains. A volume fraction /i of 10~^ associ- 
ated with the blocking domain boundaries would account 
for the large observed value of eo-"^ Raevski et al}^ have 
obtained dielectric responses similar to those of CCTO 
in the nonferroelectric perovskite ceramics AFei/2Bi/203 
with A being Ba. Sr, or Ca, and B being Nb, Ta, or Sb. 
They interpret their results as indicating that a Maxwell- 
Wagner mechanism^*^ is responsible for the large permit- 
tivity and its temperature dependence, and they use a 
simple layered structure to illustrate their argument em- 
pirically. The layered geometry is certainly too artificial, 
however, as in reality CCTO is likely to possess a more 
complicated (and strongly nonplanar) arrangement of do- 
mains and boundaries. 

A less unrealistic morphology would be spherical con- 
ducting domains surrounded by spherical shells of in- 
sulator, with e(cL') derived within the effective medium 
approximation.^^ The result is 



Sccjoj) - 2{ec{uj) - eb) fi 
3eb + {ec{<^) - eb)fi 



-£6 



where 



ec{iL)) = eb + ^mab{oj)/oJ 



(20) 



(21) 



is the complex dielectric constant of the conducting do- 
mains, and wc have set ei — ey, for simplicity. This can 
also be rewritten in the Debyc form (17) with 



eo = 3(1 - fi)eb/fi , 
t{u) = eo/(l - fi)4Trab{u) . 



(22) 
(23) 



Note that eo is increased above the lower bound in (16) 
by a factor of 3(1 — /i) in this morphology, according to 
the effective medium approximation. In both morpholo- 
gies, eo increases monotonically as /i decreases, a feature 
common to all nonpercolative morphologies. 

The above simple treatments lead to deviations from 
Debye relaxation only through the frequency dependence 
of (Jb{co), which may or may not be significant for these 
materials given our lack of understanding of the conduc- 
tion mechanisms. In particular, randomness in the mor- 
phology can lead to such deviations, particularly as close 
to the percolation threshold as the samples appear to be, 
with /i ~3xl0~^ for eb ~10^ and eo ~10^. Our careful 
examination of the observed frequency dependence of e' 
and e" shows that the expected deviations from Debyc- 
like behavior are indeed present.'' This case of conduct- 
ing domains enveloped by insulating domain boundaries 
is identical to that of cermets with a low volume fraction 
of ceramic or of varistors. Since these have been exhaus- 
tively studied, we need comment no further on this case. 
Note that Eqs. (17), (22), and (23) provide a basis for 
the generic form of ej^, Eq.(9), introduced without justi- 
fication in Sec. IIIB. 



VI. UNBLOCKED MOPHOLOGIES WITH 
CONDUCTING BOUNDARIES 

In this section, we present a model of a spatially in- 
homegenous material whose bulk is insulating but whose 
internal boundary or interface regions are conducting. 
The model is developed from the premise that the con- 
ducting boundaries are randomly parallel to one of three 
orthogonal cubic planes; as a simplification, we initially 
treat the boundaries as disjoint disks which we further 
simplify to oblate ellipsoids, identical except for their 
random orientation, with principal axes a = h ^ c and 
occupying volume fraction (f). For such a model to ap- 
ply to CCTO, continuous conducting pathways through 
the sample formed from boundaries are forbidden, and 
so, while the conducting boundaries are allowed to in- 
tersect, we force the sample to be in a nonpercolating 
regime. The complex dielectric constant e* of our sim- 
ple morphogenetic model system, derived below via an 
effective medium approximation, agrees quantitatively 
with experiments, reproducing both the magnitude and 
frequency-dependence of the observed dielectric response 
in CCTO and related materials. 

We begin by supposing that our model material con- 
tains a density N of conducting boundary segments, 
which for simplicity wc take to be in the form of dis- 
joint thin discs of radius a' and thickness 2c' randomly 
distributed in the three cubic planes with volume fraction 



27r A^a'^c'. 



(24) 



Since the overall symmetry is still cubic, we note that 
e is a scalar, and thus the response can be computed 
by assuming an electric field and polarization parallel to 
a single cubic axis. We now derive e* via an effective- 
medium approximation. Suppose that each disc is at the 
center of a sphere of radius h of the bulk material with 
complex susceptibility Xhi where 



eb = 1 -I- 47rxb, 



(25) 



outside of which is the effective medium. The volume 
fraction within the sphere of bulk material then becomes 



3 a'V 
2 63 



(26) 



which leads to an important condition: for the conduct- 
ing discs to contain the spheres completely and there- 
fore be below the percolation threshold, we must have 
< 1 in order to make contact with experiments. 



A. Bulk polarization 

The total macroscopic polarization contains contribu- 
tions from the bulk-like insulating domains, and also 
their associated internal conducting interfaces. We de- 
rive each contribution in turn. The polarization in the 
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bulk material, Ph, is related to the local field within the 



bulk, E^"", as 



(27) 



the local field inside the bulk material can be expressed 
in terms of local and macroscopic polarizations, and also 
the macroscopic field, as 



loc 

E^'^ = Ejna,c H ^(-fmac ~ Ph) t 



(28) 



where -Emac is the macroscopic field in the effective 
medium and Pmac is given by 

e* - 1 

-fmac — Xmac -^mac — Z -^mac ) (2^) 



where 



47r 



e* = 1 + 47rXmac • 



(30) 



To obtain e*, we need an expression for the macro- 
scopic polarization, Pmac, in terms of .Emac- We begin by 
first obtaining the former, Pb, in terms of Smac- Inserting 
(29) into (28) gives 



An 



(31) 



now, inserting (31) into (27), we obtain the desired result 

i(e*+2)xb 



Ph 



47r , -Emac 



1 + 



(32) 



B. Polarization of the conducting internal interfaces 

Deriving the conducting-interface contribution to the 
macroscopic polarization is slightly more involved, and, 
to simplify the discussion without loss of generality, we 
now replace conducting discs with ellipsoids of the same 
volume 47ra^c/3 = 27ra'^c' and the same eccentricity 



- 1 » 1, 



(33) 



where a = (3/2)^/'^a', so that the polarization and elec- 
trostatic field inside remain uniform. The volume frac- 
tion of ellipsoids to bulk material (Eq. (26)) then becomes 



The local field E^°^ inside the ellipsoid differs from _E^°^ 
by the depolarizing field introduced by the surface charge 
arising from the change in polarization across the the 
ellipsoid-bulk interface: 



^iT = 4°' + ia(^b-f"ai), (35) 

where P„, iw the polarization inside, L,,. is the depolar- 
ization factor, and a specifies the cubic axess normal 
to the boundary plane, i.e., normal to the semi- minor 
axis c, where La = An. Using the linear relation 
P^i = ;Yr^aT. and also Eq. (35), we find 



where 



X*(l + LaXb) 

1 + LaX* 



Xi = Xb + — 



TTiloc 



(36) 



(37) 



is the complex susceptibility of the boundary region, Ui 

its conductivity, and the real part of x* has been taken 
as Xb for simplicity. Using (31) and (32), E^^ can be 
expressed in terms of E'mac as 



^loc _3 



Ue*+2)E„ 



1 + f Xb 



(38) 



and from (36) and (38), we relate Pai and .Emac, as de- 
sired: 



P.. = l(e + 2)^*(1 + ^"^^) 



3'- ■ 1 + f Xb 



(39) 



We now impose the cfFectivc-mcdium self-consistency 
condition that Pmac is the volumetric average of Pb, i.e., 



Pmac = (1 - ^)Ph + ' 



(40) 



where we have taken the orientational average Pj = 
^J2aPai- The above equation, with the values of Pd 
and Pb derived above, expresses Pmac linearly in terms 

of E^mac- 



C. Dielectric response 

We now recast e* = 1 +47rPmac/-Emac in the suggestive 
form 



, _ Cb + ^BjCi 



where we define 



1 1 



(34) and 



■ IWTr, 



(41) 



(42) 



(43) 



The quantities Pi can be thought of as inverse depolar- 
ization factors, with an additional dynamical correction. 
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Equation (41) is a generalization of the Clausius-Mossotti 
relation which displays the proper limiting behavior, as 
follows. First, e reduces to eb when (f> vanishes. Second, 
it reduces to the Clausius-Mossotti relation in standard 
form for an array of metallic spheres for which = 0, 
Bi — > 3/47r as w ^ 0; that is, 



1 + 

1 - ^Na' 



(44) 



where a = is the polarizability for the case that c = a. 
Third, when (Ti — > 0, the material becomes uniform and 
e* reduces to eb- 

In the present case, Bj is the sum of two terms and can 
be written explicitly as 



2 1 



1 



1 1 



1 



3 La 1 — ioJTa '3 Lcl — ioJTc 



Ba + B„ (45) 



where Ba and Be are associated with the semimajor axes 
a and c, and 



(Xb + l/La)/ai 



and 



Tc = (Xb + l/Lc)/ai. 
Prom Landau and Lifshitz,^^ we obtain 



-(e — tan ^ e) 



La = -(47r-Lc) , 



(46) 



(47) 



(48) 



(49) 



with e given by (33). For the oblate ellipsoids, a/c and 
e are large enough so that Lc approaches An and La be- 
comes small, 



La=7r' 



(50) 



Inserting Lc = 47r and (50) for La into (45) and (46) 
yields 



5a = 



2 a 1 



and 



with 



and 



B, 



c 1 — iuiTa 



1 1 
127r 1 — iuiTc 



If + eb + l 



Tb 



47r(Ti 



eb 
47rcri 



(51) 



(52) 



(53) 



(54) 



Ba and Be assume their maximum magnitudes at a; = 0; 
explicitly. 



(55) 



provided a/c^ 1. Since </> <C and eb ~ 80-90 (Refs. 3, 
7), (87r/3)i?c(0)(/) is negligible relative to eb in the numer- 
ator of e*, (41), and (47r/3)-Bc(0)?!> is negligible compared 
with unity in the denominator of e*. Ftom (51), the ratio 
of the relaxation times is 



Tb TT ebC eb 



(56) 



Under the suppositions that a 3> c and <C 1, it can be 
shown from (51), (53), and (56) that Be docs not con- 
tribute significantly to e* at any frequency, independent 
of the relative magnitudes of a/c and Cb in (56). 

Be can thus be safely neglected, and e* can be put into 
the Debye form with eoo = eb and 



eo = 



eb + ^Na^ 
1 - ^Ncfi ■ 



(57) 



Remarkably, Eq. (57) allows the dielectric constant to be 
interpreted as that arising from a collection of metallic 
spheres of effective radius 



97r 



1/3 



(58) 



embedded in a medium of dielectric constant eb- 
Equation (57) can be rewritten as 



eo 



eb + 2/ 
1-/ ' 



where 



3 



(59) 



(60) 



is the volume fraction of the regions of the sample within 
which the electrostatic field is screened out by the a- 
oriented conducting boundaries, an analogue of the Fara- 
day cage effect. The Debye relaxation time is now given 
by 

(61) 



1-/ 

As for Eq. (23) for the NP/CB/UB case, aU temperature 
and frequency dependence of r arises from the conduc- 
tivity (Ti. 

There is a pseudo percolation threshold at / = 1 
at which e in Eq. (59) diverges when the field-free re- 
gions percolate. However, this mean-field theory and 
the geometry of the model both grossly overestimate 
the value of the percolation threshold through neglect 
of randomness.^® This overestimated threshold occurs 
at a value of (/> about 1.5% below the true percolation 
threshold of our disk-in-a-sphere model when the disk 
and sphere radii coincide, also overestimated. 
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VII. DISCUSSION AND CONCLUSIONS 

We have constructed six possible morphologies for in- 
homogcncity in the local dielectric response of CCTO 
by supposing the properties of the internal boundary re- 
gions to be opposite from those of the bulk. (If one 
conducts, the other is insulating, and vice versa.) We 
have pointed out that the large dielectric constants re- 
sult either when the conducting regions almost percolate 
or when they percolate but are blocked at the surface 
or electrode interface. We have analyzed the expected 
dielectric properties of all six morphologies in Sees. III- 
VI, in some cases explicitly and in others by model cal- 
culation. Without exception, the dielectric properties of 
each of the six morphologies can be consistent with all re- 
ported observations. ^^"^ Our principal conclusion is that 
existing experiments do not distinguish between grossly 
different morphologies, e.g., conducting bulk inhibited 
from percolation by insulating internal boundaries, and 
conducting internal boundaries embedded in and dis- 
rupted by insulating bulk. As discussed in Sec. Ill, dis- 
tinguishing between intrinsic and electrode-induced sur- 
face barriers can be achieved by studying the dielectric 
response as a function of the thicknesses of insulating lay- 
ers interposed between the CCTO sample and electrodes. 
However, the central questions are about the morphol- 
ogy of the conducting regions: are they bulk or are they 
associated with specific internal 2D faults? The nonlin- 
ear current-voltage characteristics of the two morpholgies 
should differ quantitatively. In both cases there should 
be a threshold field at which current rises rapidly with 
voltage, going over to a superlinear voltage dependence. 
In the conducting boundary morphology (case 3), field 
concentration would occur at the breaks in the current 
pathways along the boundaries; this would not occur if 
the bulk were conducting (case 2), and thus its threshold 
voltage would be higher. 

Additional dielectric measurements at low frequencies 
(<1 MHz) and temperatures (T<100 K) may be able 
to distinguish morphologies with conducting bulk (CB) 
from those with conducting boundaries or interfaces (CI). 
According to Homes et al.,^ the intrinsic bulk dielectric 
constant eo is enhanced by about 75% as the temperature 
is reduced to 10 K from room temperature, a surprising 
increase which can be linked mainly to an increase in os- 
cillator strength of one low-frequency IR-active mode.^ 
This temperature dependence would carry over to the ex- 
trinsic giant dielectric response in CI cases, since in these 
scenarios the giant response comes about from an ampli- 
fication of the bulk cq. In CB morphologies, however, the 
intrinisic boundary eo may or may not be different from 
that of the bulk, and would not necessarily exhibit a tem- 
perature dependence. Thus, if measurements of the giant 
response at low temperatures (and lower frequencies) did 
not show the increase expected from that of the bulk eo, 
one could infer that the CI morpliologies do not occur and 
that in the CB morphologies the boundary eo differs sig- 



nificantly in its temperature dependence from that of the 
bulk. In addition, we note in both the CCTO and Cd- 
CTO eo increases dramatically at higher temperatures. 
This would follow directly from a thermally-stimulated 
increase in conductivity locally throughout the material, 
bringing it closer to the perocolation threshold. How- 
ever, such would be the case for all morphologies and, 
accordingly, is not diagnostic. 

Additionally, what is the mechanism responsible for 
the conductivity? The observed activation energy of the 
relaxation time r is that of the conductivity. The values 
obtained are consistent with extrinsic conductivity asso- 
ciated with shallow traps. It would be of interest to in- 
ject charge at one interface and measure the distribution 
of arrival times at the other to see whether anomalous 
transport occurs. 

Observations with probes sensitive to the conductivity, 
stiTicture, and composition at the nanoscale would be of 
significant experimental interest. In Sec. II we discussed 
a possible mechanism by which conductivity might arise 
in a Cu-Cu antiphase boundary. In (Ba,Sr)Ti03 A-A and 
B-B antiphase boundaries, Naumov et al.^^ have found 
compressive stress via first-principles computations. This 
could be relaxed by a locally increased concentration of 
oxygen vacancies leading to activated n-type conduction 
in the boundaries, as discussed in Sec. I for vacancies in 
the bulk. Measurements of the sign of the thermopower 
of samples with measureable conductivities could there- 
fore be interesting. However, one should recognize that 
the conducting boundary morphologies require fine tun- 
ing of the boundary volume fraction to remain close but 
not exceed the percolation threshold. 

Finally it is important to recognize that the mecha- 
nisms leading to large dielectric constants can differ in 
polycrystalline (ceramic) and single-crystal samples. For 
the former, a model of conducting grains and blocking 
grain boundaries may be the most plausible since, if the 
reverse were true, the grain boundaries would percolate 
and there would be substantial conductivity. For the lat- 
ter, we cannot yet distinguish between conducting and 
insulating internal boundaries. Tselev et al.,^'^ however, 
have argued for the conducting bulk and insulating inter- 
nal boundary morphology, the conductivity being associ- 
ated with oxygen vacancies and the insulating boundaries 
being deficient in oxygen vacancies. 
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